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ABSTRACT 

Accurate  performance  evaluation  has  always  been  an 
important  issue  in  network  design  and  analysis.  Dis¬ 
crete  event  simulation  has  been  known  to  be  accu¬ 
rate  but  very  time  consuming.  Thus  analytical  meth¬ 
ods/approximation  is  necessary  for  real  time  estimation, 
large  scale  network  optimization  and  sensitivity  analy¬ 
sis.  In  a  circuit-switched  loss  network,  a  particular  per¬ 
formance  metric  of  interest  is  the  end-to-end  blocking 
probability.  Various  analytical  approaches  and  approx¬ 
imation  schemes  have  been  suggested  for  this  problem 
and  among  them.,  the  fixed-point  method,  or  reduced  load 
method,  has  received  much  attention.  However,  most  of 
these  schemes  considered  either  only  single  traffic  rate 
situation  or  multi-rate  traffic  under  fixed  routing.  We 
develop  an  approximation  scheme  to  estimate  end-to- 
end  blocking  probability  in  a  multi-rate  multi-hop  net¬ 
work  with  an  adaptive  routing  scheme.  The  approxima¬ 
tion  results  are  compared  with  that  of  discrete  event  sim¬ 
ulation.  An  example  of  application  is  also  provided  in 
which  the  proposed  scheme  is  linked  to  the  optimization 
tool  CONSOL-OPTCAD  to  get  network  design  trade¬ 
offs.  The  method  described  here  is  readily  applicable 
to  the  accurate  performance  evaluation  of  military  net¬ 
works  which  are  often  large  and  hybrid. 

INTRODUCTION 

This  paper  is  focused  on  the  evaluation  of  one  network 
performance  metric  -  the  call  blocking  probability  of  a 
loss  network  [1], 

There  are  two  types  of  approaches  to  evaluate  the  call 
blocking  probabilities,  or  to  evaluate  other  performance 

*  Prepared  through  collaborative  participation  in  the  Advanced 
Telecommunications /Information  Distribution  Research  Program 
(ATIRP)  Consortium  sponsored  by  the  U.S.  Army  Research  Lab¬ 
oratory  under  the  Federal  Research  Laboratory  Program,  cooper¬ 
ative  agreement  DAAL01-96-0002  and  under  NASA  cooperative 
agreement  NCC3-528. 


metrics.  One  is  some  sort  of  analytical  approach  (ap¬ 
proximation,  estimation)  and  the  other  is  discrete  event 
simulation.  Discrete  event  simulation  is  a  widely  used 
and  very  helpful  tool  of  accurate  performance  evalua¬ 
tion.  However,  even  for  a  fairly  small  network  model, 
it  takes  considerable  amount  of  run  time  to  get  satis¬ 
factory  results.  With  discrete  event  simulation,  it  is 
very  hard  to  get  the  sensitivity  analysis  of  the  perfor¬ 
mance  and  optimization  schemes,  for  which  numerous 
simulation  runs  are  needed.  Military  networks  are  typi¬ 
cal  hybrid  networks,  containing  satellites  and  terrestrial 
wireless  subnetworks.  These  networks  are  also  typically 
large.  Quite  often  we  need  to  estimate  the  performance 
of  a  network  prior  to  acquisition.  This  cannot  be  done 
in  a  reasonable  amount  of  time  by  discrete  event  simula¬ 
tion.  Another  common  problem  is  to  estimate  the  per¬ 
formance  of  a  small  network  within  a  larger  network. 
One  can  do  a  discrete  event  simulation  of  the  smaller 
network,  while  some  analytical  way  is  needed  to  repre¬ 
sent  the  effects  of  the  coupling  and  of  the  larger  network 
in  an  aggregate  fashion.  For  instance,  the  small  network 
can  be  two  wireless  LANs  serving  two  platoons  of  sol¬ 
diers,  while  the  larger  network  may  involve  many  more 
nodes  and  satellites  and  UAVs. 

Therefore  this  makes  analytical  approaches  essential, 
since  they  are  generally  much  faster  than  their  simu¬ 
lation  counterpart.  In  1917  the  Danish  mathematician 
A.  K.  Erlang  published  his  famous  formula  [2]  which  es¬ 
timates  the  loss  probability  of  a  conventional  telephone 
network.  It  covers  the  case  of  a  single  link  with  calls 
at  a  single  rate.  Analytically,  when  there  are  multi¬ 
ple  links  and  multiple  call  rates,  with  different  band¬ 
width  requirements  and  with  a  fixed  route  associated 
with  a  certain  source-destination  pair,  a  loss  network 
is  modeled  as  a  multidimensional  Markov  process  with 
the  dimension  of  the  state  space  of  the  process  being 
the  product  of  the  number  of  routes  allowed  in  the  net¬ 
work  and  the  number  of  service  classes  with  different 
bandwidth  requirements.  When  alternative  routes  are 
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present  in  addition  to  fixed  routes,  the  Markov  process 
no  longer  bears  a  product  form,  and  the  equilibrium 
state  probabilities  can  be  obtained  by  writing  out  the 
whole  set  of  detailed  balance  equations  and  solving  it  [2], 
However,  this  approach  is  not  practical  in  dealing  with 
large  networks  with  hundreds  of  thousands  of  routes 
and  integrated  services  with  multiple  service  rates,  since 
the  computational  complexity  is  both  exponential  in  the 
number  of  routes  and  exponential  in  the  number  of  ser¬ 
vice  classes.  This  leads  to  the  need  for  the  development 
of  computational  techniques  that  provide  accurate  esti¬ 
mates  of  blocking  probabilities  for  loss  networks. 

The  reduced  load  approximation,  also  called  the  Erlang 
fixed-point  method,  has  been  proposed  for  this  scenario 
and  has  received  much  attention  [2],  The  reduced  load 
approximation  is  based  on  two  assumptions:  One  is  the 
link  independence  assumption  which  assumes  that 
the  blocking  occurs  independently  from  link  to  link.  The 
other  assumption  is  the  Poisson  assumption  which 
assumes  that  the  traffic  flow  to  each  individual  link  is 
Poisson  and  that  the  corresponding  traffic  rate  is  the 
original  external  offered  rate  thinned  by  blocking  on 
other  links  of  the  path,  thus  called  the  reduced  load. 

Most  of  the  earlier  works  in  fixed-point  method  either 
studied  the  multi-rate  situation  with  fixed  routing  [5], 
or  focused  on  state-dependent  routing  schemes  with  sin¬ 
gle  traffic  rate  [5],  or  multi-rate  service  with  single  link 
(resource)  [6].  In  [4],  a  fixed-point  method  was  proposed 
to  approximate  blocking  probabilities  in  a  multi-rate 
multi-hop  network  but  introduced  additional  computa¬ 
tional  efforts  in  solving  the  associated  network  reliability 
problem  are  needed,  and  also,  the  routing  scheme  used 
was  sequential  routing  with  trunk  reservation,  which  is 
not  often  practiced  in  real  networks. 

We  focus  our  attention  on  the  evolving  integrated  ser¬ 
vice  networks  which  bear  the  following  characteristics: 
(1)  The  networks  are  typically  much  sparser  and  have  a 
more  hierarchical  topology.  Thus,  the  assumption  of  the 
existence  of  a  direct  route  between  source  and  destina¬ 
tion  nodes  does  not  hold  in  most  instances.  (2)  Routes 
can  comprise  a  much  larger  number  of  hops  (typically 
around  5  or  6)  and  there  are  typically  a  large  number 
of  possible  routes  between  source  and  destination  nodes. 
(3)  The  presence  of  different  traffic  classes  characterized 
by  widely  varying  bandwidth  requirements  and  different 
mean  holding  times  must  be  considered. 

Motivated  by  the  above,  we  propose  to  use  adaptive 
routing  in  combination  with  the  fixed-point  method  to 
calculate  call  blocking  probabilities. 


A  FIXED-POINT  METHOD  FOR 
BLOCKING  PROBABILITY  ESTIMATION 

Consider  a  network  with  J  links,  each  indexed  as  ('«,  j), 
with  i .  j  being  the  index  of  the  link  end  nodes.  Link 
(z,  j )  has  a  capacity  of  C\:1  units  of  circuits  or  trunks. 
The  network  supports  a  total  of  S  classes  of  connections, 
where  a  class-5  connection  has  a  bandwidth  requirement 
of  bs  G  on  every  link  on  the  path  the  connection  is 
routed.  Every  class-5  connection  also  has  an  arrival  rate 
of  A(s),  and  the  mean  time  of  call  duration  is  l/n(s). 

For  each  source-destination  node  pair  ( k,l ),  there  is  an 
associated  set  of  routes  7 Zkj  =  ...,TZpp}^ 

with  each  route  being  a  subset  of  the  set  of  links: 

£  {1,...,J}.  A  connection  is  accepted  if  some 
route  has  available  bandwidth  on  each  of  its  links  to 
accommodate  this  connection,  and  the  connection  is 
routed  on  that  route  and  holds  the  bandwidth  for  a 
duration  with  mean  time  1/7/ (5).  If  none  of  the  routes 
are  available,  the  connection  is  rejected.  The  end-to-end 
blocking  probability  of  a  class-5  connection  with  source 
destination  pair  ( k,l )  is  denoted  by  Bkj(s).  Through¬ 
out  this  paper  the  link  is  considered  to  be  duplex  and 
bi-directional. 

The  common  routing  policies  which  have  been  studied 
are  fixed  routing,  alternative  routing,  sequential  alter¬ 
native  routing  and  adaptive  alternative  routing.  We 
focus  on  the  last.  One  important  scheme  of  this  kind  is 
called  the  Least  Loaded  Routing  (LLR),  where  the 
call  is  first  tried  on  the  direct  route,  if  there  is  one.  If 
it  cannot  be  setup  along  the  direct  route,  the  two-link 
alternative  route  with  the  largest  number  of  point-to- 
point  free  circuits  is  chosen.  A  version  of  LLR  was 
recently  implemented  in  the  AT&T  long-distance  do¬ 
mestic  network  [2],  The  method  we  investigated  here  is 
an  extension  of  LLR  and  is  also  motivated  by  [3]  and 
[4].  It  is  a  min-max  scheme  in  that  it  tries  to  maximize 
the  free  bandwidth  on  the  link  which  has  the  minimum 
free  bandwidth  on  a  route.  Call  admission  control  of 
the  trunk  reservation  type  is  also  considered. 

Maximal  Residual  Capacity  Adaptive  Routing 
Scheme 

The  routing  policy  considered  here  is  described  as  fol¬ 
lows.  Each  source-destination  node  pair  ( k ,  l)  is  given  a 
list  of  alternative  routes  =  {IZpf},  m  =  1,2,...,  M. 
When  a  call  arrives,  each  of  the  routes  on  the  list  is 
evaluated  to  determine  the  number  of  free  circuits  on 
each  link  employed  by  the  routes. 

Let  Op/6  denote  the  free/available  bandwidth  on  link 
(i,j)  when  the  call  arrives,  and  consider  a  route  TZpf 
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which  comprises  of  N  links  expressed  by  = 

{(*1,  ji),  (*2,  J2),  •••,  (*w,  Jn)}-  The  route  is  in  a  state 
of  admitting  a  call  of  bandwidth  requirement  bs  only  if 
C/’f  >  be,  n  =  1,2 . N  under  no  trunk  reservation 

admission  control  or  C/,ef  >6,  +  rs,  n  =  1,2, ....  N 
under  trunk  reservation  admission  control,  where  rs  is 
the  trunk  reservation  parameter  of  the  call. 

Consider  a  route  Tiffi  which  is  presently  available  for 
a  call,  the  most  congested  link  on  the  route  is 

the  link  with  the  fewest  free  circuits  on  this  route: 
44  =  arqmin,.  .  ,  ^tm)Cfieie .  When  there  are  more 

('hi,Jn)enkl  ‘n,Jn 

than  one  route  available  in  the  route  list  TZk.l-  the  one 
with  the  maximum  free  bandwidth  on  its  most  con¬ 
gested  link  is  selected  for  connecting  the  call.  In  the  pro¬ 
posed  approximation,  we  consider  steady  state, 

^ k,i 

must  be  replaced  by  E[C^(^],  44  becomes  the  statis- 

^k,i 

tically  most  congested  link  as: 

=  argmax^^Zij 

where  zy  is  defined  as  the  link  load  of  link  zy  = 

J2S  4(1)  'chi  where  Vi,j{s)  is  the  reduced  load/arrival 
rate  for  class-5  calls  on  link  If  none  of  the  routes 

are  admissible,  the  call  is  blocked. 

Algorithm 

In  the  proposed  method,  the  fixed  point  is  achieved  by 
mappings  between  the  following  four  sets  of  unknown 
variable:  ;'y  (5),  ay  (5),  Pij{n)  and  q^m\k,l,s),  as  il¬ 
lustrated  in  the  figure  below: 


Figure  1:  Relationship  between  variables 

where  ay  (5)  is  the  probability  that  link  i.  j  is  in  a  state 
of  admitting  class-5  calls;  py(n)  is  the  link  occupancy 
probabilities  for  link  i.e.,  the  probability  that  ex¬ 

actly  n  units  of  circuits  are  being  occupied  on  link  (*,  j); 
and  q(m\k,  l,  s)  is  the  probability  that  a  call  of  class-5 
with  source-destination  pair  of  (k,  l)  is  attempted  on  its 
rn'h  route. 

Mapping:  aij(s),  q^m\k,l,  s)  =>■  Vij(s)  The  link  ar¬ 
rival  rates  on  link  (i,j)  from  the  node  pair  (k,l)  for 
class-5  traffic  on  the  rnlh  route  is  given  by  the  reduced 
load  approximation  as 

4?  (*>  s)  =  xk,l(s)q{m)  {k,  l ,  s)I[{i,  j)  G  nff] 


aU.V  (5  ) 

(u,v)  en^  ,(u,v)^(i,j) 

where  I  is  the  indicator  function.  The  aggregate  arrival 
rate  of  calls  of  class  5  on  link  (i,j)  is  given  by 

viAs)  =  EEi/!?^i-s) 

k,l  m 


Mapping:  vitj{s)  =k  ay(5),py(n) 

Given  ty  (5),  we  can  compute  the  link  occupancy  prob¬ 
abilities  Pij(n)  for  each  link  in  the  network.  This  can 
be  done  by  either  using  Kaufman’s  simple  recursion  [7] 
when  there  is  no  trunk  reservation  present,  or  using  ap¬ 
proaches  proposed  by  [6,  8]  as  suggested  by  Greenberg 
in  [4], 

Mapping:  py(n)  =^>  q(m)  {k,l,s) 

Given  py(n),  define  for  link  (i,  j),  the  probability  of  no 
more  than  n  trunks  are  free  (at  most  n  trunks  are  free) 
as:  7y(n)  =  Ylt=oPi,j^i,j  ~  ^)-  So  the  probability  of 
attempting  the  call  on  the  rnlh  route  is: 


c 


(m) 


qim)(k,l,s)  = 


k,l 

Up, 

n= 0 


m— 1 

(m)  (C  (m)  —  fl)  JJ  T^(t){n  ~  1) 

k,l  ^ k,l  _ j  ^ k,l 


M 


II  T,«(n) 


t=m+ 1 


k,l 


End-to-end  Blocking  Probabilities 

Finally,  the  end-to-end  blocking  probability  for  calls  of 
class  5  between  source-destination  node  pair  (k,l)  is 
given  by 

B{k,l,s)  =  1  -J2q(m)(k,l,s)  ««,■»(■§)• 

m  Men™ 

Repeated  substitution  should  be  used  to  obtain  the 
equilibrium  fixed  point.  And  finally  the  end-to-end 
blocking  probabilities  can  be  calculated  from  the  fixed 
point. 


EXPERIMENTS  AND  EVALUATION 


This  example  is  borrowed  from  [4]  with  minor  changes. 
The  topology  is  derived  from  an  existing  commercial 
network  and  is  depicted  in  Figure  2  below.  The  discrete 
event  simulation  was  done  by  using  OPtimized  Network 
Engineering  Tools  (OPNET). 

There  are  16  nodes  and  31  links,  with  link  capacity 
ranging  from  60  to  180  trunks.  The  detailed  link-by¬ 
link  traffic  statistics  and  link  capacities  can  be  found 
in  either  [4]  or  [10].  The  traffic  in  the  network  consists 
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Figure  2:  Topology  of  Example  Network 


of  four  types,  namely  class- 1,  2,  3,  and  4,  and  require 
bandwidth  of  1,  2,  3,  and  4  trunks,  respectively.  No 
admission  control  is  employed. 


Node  Pair 

Class 

FPA 

DES 

(0,4) 

4 

0.018213 

(0.0122,  0.0179) 

(0,13) 

1 

0.074434 

(0.0729,  0.0766) 

(1,6) 

1 

0.077371 

(0.0697,  0.0701) 

(5,6) 

3 

0.229528 

(0.2262,  0.2278) 

(6,10) 

2 

0.147436 

(0.1420,  0.1483) 

(9,13) 

4 

0.307191 

(0.2794,  0.2848) 

Number  of  Iterations 

28 

CPU  Time(seconds) 

145.43 

4.3  x  104 

Table  2: 

1.4  Times  The  Nominal  Traffic. 

Node  Pair 

Class 

FPA 

DES 

(0,4) 

4 

0.112658 

(0.0025,  0.0026) 

(0,13) 

1 

0.135564 

(0.1492,  0.1500) 

(1,6) 

1 

0.156322 

(0.1445,  0.1466) 

(5,6) 

3 

0.419781 

(0.3922,  0.3940) 

(6,10) 

2 

0.269145 

(0.2572,  0.2583) 

(9,13) 

4 

0.519083 

(0.4791,  0.4793) 

Number  of  Iterations 

24 

CPU  Time(seconds) 

125.11 

2.3  x  106 

Table  3:  1.8  Times  The  Nominal  Traffic. 


APPLICATION  IN  NETWORK  DESIGN 
USING  CONSOL-OPTCAD 


In  routing,  any  node  pair  is  allowed  routes  that  have  at 
most  4  hops.  Multiple  routes  for  one  node  pair  are  listed 
in  order  of  increasing  hops,  with  ties  broken  randomly. 
Each  link  is  considered  to  have  same  unit  length,  so  only 
the  hop  number  is  counted. 

Results  for  some  selected  node  pairs  and  classes  are 
listed  in  Table  1  through  3(FPA  stands  for  fixed-point 
approximation  and  DES  stands  for  discrete  event  sim¬ 
ulation),  each  corresponding  to  a  different  traffic  load. 
Table  1  corresponds  to  the  “nominal”  traffic  which  is 
provided  in  [4]  and  [10].  Tables  2,  3  show  the  results  for 
traffic  1.4  and  1.8  times  the  nominal  traffic,  respectively. 


The  proposed  fixed-point  approximation  gives  conserva¬ 
tive  estimations  generally,  and  it  improves  as  the  load 
gets  heavier.  These  results  strengthen  the  argument 
that  these  approximations  are  indeed  very  useful  as  es¬ 
timators  of  worst  case  performance. 


Node  Pair 

Class 

FPA 

DES 

(0,4) 

4 

0.000178 

(0.0,  0.0) 

(0,13) 

1 

0.006341 

(0.0021,  0.0034) 

(1,6) 

1 

0.006473 

(0.0030,  0.0034) 

(5,6) 

3 

0.020463 

(0.0189,  0.0201) 

(6,10) 

2 

0.013222 

(0.0109,  0.0138) 

(9,13) 

4 

0.028468 

(0.0185,  0.0245) 

Number  of  Iterations 

18 

CPU  Time(seconds) 

94.1 

3.7  x  104 

Table  1:  Nominal  Traffic. 


One  of  the  main  reasons  we  are  interested  in  devel¬ 
oping  an  analytical  approximation  algorithm  is  be¬ 
cause  such  an  algorithm  can  be  easily  linked  to  math¬ 
ematical  programming  tools  to  get  network  perfor¬ 
mance  optimization  and  trade-off  analysis.  In  this  sec¬ 
tion,  we  link  the  proposed  reduced  load  approximation 
method  with  CONSOL-OPTCAD  [9],  which  is  a  tool 
for  optimization-based  design,  and  show  how  network 
parameter  design  can  be  realized. 

The  application  we  have  here  is  to  design  the  trunk 
reservation  parameters.  As  a  way  of  call  admission  con¬ 
trol,  trunk  reservation  regulates  individual  classes  of 
traffic  as  well  as  their  interrelationship.  How  to  choose 
the  combination  of  {ri,  r2, ...,  r,?},  where  rs  is  the  trunk 
reservation  parameter  of  class-5  traffic  and  S  is  the  total 
number  of  different  classes  the  network  carries,  to  get 
the  best  network  performance  is  important.  Here  the 
network  performance  of  interest  is  the  average  blocking 
probability  of  the  network  as  a  whole  and  the  blocking 
probabilities  that  each  individual  class  of  traffic  experi¬ 
ences. 

This  design  problem  is  formulated  as  follows:  Design 
parameters:  r\,  r27 ...,  rs-  Objective: 


Minimize 


S(fc,Q,s  ^kAs)  '  [1  ~  Bk,l(s)]  '  Bk,lis) 
E(fe,i),s  A kj.(s)  •  [1  -  Bkj{s)] 
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Constraints: 

Bk,i(s)  <  boundk,i{s),  all  ( k,l )  and  s 

where  B}.,i(s)  is  the  blocking  probability  of  class-s  traf¬ 
fic  between  source-destination  node  pair  (k,l),  and 
boundkj(s)  is  the  upper  bound  for  this  probability. 

By  applying  this  model  to  a  fully  connected  five-node 
network  serving  3  classes  of  traffic  [10]  with  trunk  reser¬ 
vation  admission  control,  we  get  the  trade-off  between 
the  blocking  probability  of  each  class  vs.  the  weighted 
average  blocking  probability,  which  is  shown  in  Table  4 
and  Figure  3  below  (The  uppermost,  middle  and  bot¬ 
tom  curves  are  blocking  probabilities  of  Class  3,  2  and 
1,  respectively). 

As  we  see  from  Table  4,  the  weighted  average  block¬ 
ing  probability  and  the  blocking  probability  of  class- 1 
type  of  traffic  achieve  their  optimum  at  the  same  time 
with  trunk  reservation  parameter  choice  of  1,4  and  5. 
The  reason  is  obvious:  since  class- 1  has  much  smaller 
trunk  reservation  requirement  than  class-2  and  3,  to¬ 
gether  with  its  lowest  bandwidth  requirement,  it  has 
the  highest  priority  and  chances  of  being  admitted  into 
the  network.  On  the  other  hand,  class-2  and  3  are  be¬ 
ing  jeopardized  by  their  high  trunk  reservation  require¬ 
ment  and  also  high  bandwidth  requirement.  This  phe¬ 
nomenon  can  also  be  observed  in  Figure  3.  The  curves 
seem  very  random  here  because  they  are  basically  the 
connection  of  discrete  samples  at  different  trunk  reser¬ 
vation  parameter  combination  choices. 


Average 

Bi 

b2 

b3 

(ri,r2,r3) 

0.266 

0.036 

0.580 

0.803 

(1,4,5) 

0.273 

0.063 

0.527 

0.841 

(1,3,5) 

0.290 

0.068 

0.563 

0.792 

(2,4,5) 

0.293 

0.112 

0.467 

0.861 

(1,2,5) 

0.307 

0.117 

0.486 

0.819 

(1,2,4) 

0.310 

0.120 

0.492 

0.824 

(2,3,5) 

0.420 

0.571 

0.284 

0.565 

(4,1,2) 

0.406 

0.319 

0.672 

0.323 

(3,5,1) 

Table  4:  Trunk  Reservation  Parameter  Design 


CONCLUSION1 

In  this  paper  we  presented  an  approximation  scheme  of 
calculating  the  end-to-end,  class-by-class  blocking  prob¬ 
ability  of  a  loss  network  with  multi-rate  traffic  and  adap¬ 
tive  routing  scheme.  It  provides  fairly  good  estimates  of 
call  blocking  probabilities  under  normal  and  heavy  traf¬ 
fic.  We  also  presented  applications,  where  the  proposed 

lrThe  views  and  the  conclusions  expressed  in  this  paper  are 
those  of  the  authors  and  should  not  be  interpreted  as  represent¬ 
ing  the  official  policies,  either  expressed  or  implied  of  the  Army 
Research  Laboratory  or  the  U.S.  Government. 


scheme  is  linked  to  the  system  design  optimization  tool 
CONSOL-OPTCAD  to  get  network  parameter  trade-off 
analysis. 


Figure  3:  Trade-off  in  trunk  reservation  parameter  de¬ 
sign. 
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